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SUMMARY 

An analysis of a radioisotope electric cell that uses the beta emitter, 
cerium 144, was undertaken. The analysis, which relies on formats set forth for 
fission electric cells, takes into account the pertinent physical constraints in 
a more accurate manner than has been done previously. Results from the analysis 
are presented in tables and as representative plots. Comparisons between accu- 
rate and approximate solutions are also included. 


INTRODUCTION 

Presently there exists a growing interest in the production of large quan- 
tities of electric power at low specific weights for use in electric propulsion 
systems. Recently an electric powerplant design in the form of a battery that 
converts the energy of decaying radioisotopes directly into electric power has 
been proposed (ref. l). The most desirable feature of such a cell is its low- 
weight to power-output ratio (lb/kw). A further characteristic of the cell is 
its requirement of a colloidal-particle thrustor because of the high voltage 
generated. This thrustor, although presently in an early research stage, shows 
high theoretical efficiencies over a large range of specific impulse (ref. 2). 

In the original paper (ref. l) on space radioisotope cell designs, calcu- 
lations were performed from information available on fission electric cells 
(refs. 3 and 4). This previous study included calculations of operating volt- 
age, cell efficiency, weight per unit power output, and component temperature. 
Since the basic analysis on fission electric cells yielded a limited amount of 
pertinent information, the study was approximate and incomplete, and conse- 
quently a more exact analysis of the cell was undertaken. This analysis is con- 
cerned with a steady- state operation of radioisotope cells that use the beta- 
emitting radioisotope, cerium 144. A discussion relevant to the choice of this 
beta-emitting radioisotope appears in reference 5. The present report describes 
in detail the analytical method and presents the basic results of the computa- 
tions performed numerically on the IBM 7090 at the Lewis Research Center. Por- 
tions of these results were previously summarized in reference 6. 



ANALYSIS 


The operating principles of the radioisotope electrogenerator are as fol- 
lows: Beta particles (electrons), as depicted in sketch (a), are emitted in all 
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directions from the radioisotope fuel layer. The energetic beta particles travel 
through the fuel layer and out through a covering support layer. At this point 
beta particles with sufficient energy travel across the vacuum gap to the col- 
lector foil. In this action the emitter becomes positively charged and the col- 
lector becomes negatively charged, thus producing a potential field across the 
cell. It is in the crossing of this potential field that the kinetic energy of 
the beta particle is converted directly into electrical energy. In steady-state 
operation, the flow of current from the collector is compensated for by the in- 
ternal flow of beta particles to the collector. Further descriptions of the pro- 
posed operation of the radioisotope electrogenerator in an electrostatic propul- 
sion system are given in references 1, 5, and 6. 

Operation of the cell depends on the emission of highly energetic beta par- 
ticles by the radioisotope. In this radioactive decay process, beta particles 
are emitted over a continuous spectrum of energy. Sketch (b) shows a typical 
beta spectrum with n(E) representing the number of emitted beta particles per 
unit energy and E the energy of these beta particles. (All symbols are defined 
in the appendix. ) An outstanding feature of the beta spectrum is that it has a 
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maximum energy called the end-point energy (denoted as 
process is characterized by a unique end-point energy, 
of the beta-emission process yields the relation (ref. 
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m_C + E 

= _o 

E 


where 

p 

m 0 C rest energy of electron 

E energy, Mev 

e 0 value of e for E = E 0 

which holds for an allowed nuclear transition. This report considers only the 
2.98-Mer beta particle emitted from the radioisotope chain, cerium - praseodym- 
ium 144. It has been experimentally shown (ref. 8) that this beta emission can 
be represented by a Kurle plot and thus should follow the relation (eq. (l)) for 
an allowed nuclear transition. 

Since the rest energy of the beta particle is 0.511 Mev, well within the 
limits of the beta decaying energies considered in the analysis, the relativistic 
correction must be considered. Therefore, to attain an accurate description of 
the cell parameters, the relativistic formulations of the equations of motion for 
particles are used throughout this report. 


In addition to these physical limits, the problem of high-energy electrons 
traveling through materials must be considered. It will be assumed that the 
energetic electrons travel in a straight path (negligible scattering) and that 
they lose a fixed energy for a given travel distance (negligible straggling). 

In the present problem, these assumptions do not appear to be critical, since it 
can be argued that electron scattering and straggling will average out in the in- 
tended integrations. 

The Katz-Penfold relation (ref. 9) 


= 412R1* 265-0. 0954 In E mg/cm 2 (2) 

predicts the practical range FL (effective straight-line distance in which an 
electron loses its entire initial energy) for an electron of a fixed initial 
energy E (below 3 Mev). This is the best range-energy relation available for 
electrons and is a result of experiments (ref,. 9) of electron energy loss in 

aluminum. By multiplying by the factor (?r gyj (ref. 10) the Katz-Penfold rela- 
tion becomes applicable to elements other than aluminum. This form of the rela- 
tion is used in the present report to predict electron energy loss in the support 
and fuel layers of the radioisotope cell. 


It has been observed as a result of the present investigation that the use 
of nonrelativist ic equations in place of their relativistic counterparts and the 
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use of a linear approximation to the Katz-Penfold relation yielded results that 
compared closely with the nonapproximated results (see comparative section 
Effects of Changes in Physical Relations)* In view of this fact., the analysis 
that follows includes these approximate forms of the physical constraints. 

In the analysis of the direct nuclear electrogenerator cell, the following 
quantities are of interest: 

n number of beta particles reaching collector per second per unit emit- 

ter area 

electric efficiency of cell 

E c total energy of beta particles reaching collector per second per unit 

emitter area 

D(t,E d c ) distribution of beta particles reaching collector per second per unit 
emitter area as a function of their arrival angle and energy 

These parameters are important to the design of the electrogenerator cell. From 
plots of n against cell voltage, an operating point can be deduced. The pro- 
cedure for deriving this operating voltage is outlined in the section Number of 
Particles Reaching Collector. The importance of the electric efficiency is that 
it yields the value of the weight per unit power output for a particular cell - 
a measure of the applicability of the cell to space electric propulsion. Since 
the cell has three major parts (discussed subsequently), the emitter, the col- 
lector, and the vacuum gap, a knowledge of q and E e is all that is necessary 
in obtaining the entire distribution of energy in an operating cell. Finally, 
the quantity D(r, Ep^ c ) is important in considering the magnitude of backseat - 

tered electrons at the collector surface of the cell. 

Two cell geometries (see sketch (c)), the spherical concentric and the co- 
axial cylindrical, are considered. The cells consist of two foils - an emitter 
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foil as the inside element and a collector foil as the outside surface. The 
emitter foil contains a fuel layer of the beta-emitting radioisotope, cerium 144, 
spread uniformly over the Inside of a support foil. The collector is simply a 
single foil that encloses the electric field. 

For the purpose of the analysis, the following additional assumptions are 

made: 

(1) The curvature of the fuel layer is disregarded (very thin fuel layers, 
relative to the radii, will be considered). 

(2) The cell is operating at steady state. 

(3) End effects in the cylindrical configuration are not considered. 

(4) Secondary electron and backs cat tering effects are neglected. 

(5) Additional design structures (insulators, separators, connectors, etc. ) 
are neglected. 

All these conditions, it is hoped, can be approached in the actual cell design. 


SPHERICAL CONFIGURATION 
Number of Particles Reaching Collector 
A section of the fuel layer is depicted in sketch (d). Beta particles are 
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emitted from an arbitrary volume element dv, and a portion of them will leave 
the fuel layer through a unit segment of the surface dS. Since the charged 
beta particles are emitted uniformly in all spatial directions, the total number 
dn^ per second passing through dS and reaching the collector is 

dn’ « AS GOS 9 N dv for r < r "'j 
e 5 - max 

4rtr 1 

dn e = 0 for r > r Tnny 

where N g is the number of emitted beta particles per second per -unit volume 
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with an energy Eg, and r max (a function of the angle d) is some calculated 
distance to he evaluated later. 

The felement of volume in spherical coordinates is 

dV = r 2 sin 0 dr d0 dcp (4) 

and thus equation (3) expands to 

H 

dn' = — cos 9 sin 9 dS dr d0 dcp (5) 

For the spherical configuration, the coordinate system is best oriented as shown 
in sketch (e). 



Hence, the number of particles per second leaving the fuel layer per unit area 
and reaching the collector is 



where the limits r max and 9 mfly are the values of r and 9 beyond which 
there are no further contributions to the integral. (The meaning and the evalu- 
ation of the limits will be discussed more fully in the following paragraphs. ) 
Moreover, since the geometry being considered has spherical symmetry, equa- 
tion (6) reduces to 


n 


e 



dr sin 9 cos 9 


(?) 


Since the beta particles lose energy in crossing the electric field of the 
cell and also lose energy in traveling through the emitter foil, not all the beta 
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particles have sufficient initial energy to reach the collector. With this in 
mind, it can he seen that the limits of integration, r max ( 9 ) and 0 max: , are the 
values of r and 9 beyond which the initial energy of the beta particles is 
insufficient to carry them to the collector. 

The problem of finding r max ( 9 ) and 0 max may now be considered. Sketch 
(f) represents the electric field through which beta particles travel from emit' 
ter to collector. 


Collector 




(f) 



On introducing the following definitions, 

U kinetic energy of beta particle 

$ cell voltage 

e electronic charge 

P-j. linear momentum, tangential direction 

RpR£ inner and outer radii, respectively 

the equations of motion for the beta particle in the radial electric field of 
the cell can be written as 

Ui = Ug + e$ (conservation of energy) (8) 

RlRj. ^ = Rg^t, 2 (conservation of angular momentum) (9) 

where the subscripts 1 and 2 refer to the initial (emitter) and final (col- 
lector) values. 

To find the minimum initial kinetic energy U-| for which the beta 

particle can still reach the collector, the radial velocity at the collector 
v r g is set equal to zero. With this stipulation (v r 2 - 0) the relativistic 
form of equations (8) and (9) is; 
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where vq m; j_ n is the minimum initial velocity of the "beta particle and m D its 
mass. Eliminating Vj. (tangential velocity at the collector) yields 



is the minimum initial kinetic energy that a particle must have to cross the 
electric field for a given angle 9. 


In the nonrelativist ic formulation of the problem, equations (10a) and (lla) 
have the following counterparts: 


(if V? . ^ = i m Ti 
1 2 o 1, miny 2 o 


(10b) 


m R.v, . sin 9 = m R 0 v, 
o 1 l,min o 2 t 


(lib) 


which, on eliminating v^ between them, reduce to 


Equation (12b) corresponds to equation (12a) with ^ m Q v? ^ min j equivalent to the 

quantity Uq m q n and is more readily solvable than equation (12a) for the de- 
sired minimum beta-particle energy term. 

Row to complete the problem of obtaining the limits of integration, r TOBy ( 9 ) 
and 0 max , it is necessary only to consider the effect of energy loss in the 
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emitter foil (support and fuel layer). Sketch (g) shows the path of the beta 
particle through the emitter foils with values of energy placed at the base of 



each vector. Since the beta particle loses energy in passing a distance 
t g/ cos e through the support layer, its entrance energy Ep must be larger 
than U-j_ min by just the amount of this loss. Expressing the range-energy re- 
lation as R = f(E) means the condition (see sketch (h)) 

» f( El ) - f(0 ljmln ) (13) 

must be satisfied and can be evaluated to obtain the quantity E^. If the beta 
particle having the energy E at birth needs only the energy E^ (to reach 

H 
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the eolleetor) on entering the support layer, its maximum travel distance in the 
fuel layer r max (9) is 


W 8 > - f <V ■ f(E i ) (14) 

The limit 0 is simply the angle for which r max ( 6 ) first becomes zero. This 
procedure, including the evaluation of min (0), may be carried out separately 
for each angle 0. ’ 


The limits of integration r max ( 0 ) and 0 max derived in the preceding man- 
ner must be modified by the geometry of the fuel layer. The integration over the 
coordinate r must not exceed the boundary of the fuel layer (see sketch (i)). 



Thus, r „(0) is taken as the smaller of t»/cos 0 or r„ OY (0) (calculated) 

7 max i max 

■where t^ is the fuel- layer thickness. Because of this duality, a cutoff angle 
0^ exists that satisfies the relation 


cos 0^ 




1 < 


t 


f 


r 

max 


( 0 ) 


for 0.|. > 0 


► 


for 0^ = 0 


( 16 ) 


The maximum value of 0 must not exceed rc/2; thus, 0 max is taken as the 
smaller of h/ 2 or e max (calculated). Now the integral for n £ takes the form 
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sin 9 cos 9 dr 


(17) 
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which reduces to 


n 


€ 




r (0)sin 0 cos + t^(l - cos 0 + .)d0 
max 1 v 


where 0„„ must not exceed ir/2. 
max ' 


( 18 ) 


Modifications to Parameter Analysis 

It is appropriate at this point to consider the effects of the beta spec- 
trum on the parameter calculations. As discussed previously, beta particles are 
emitted over a range or spectrum of energies. The spectrum is continuous, rang- 
ing from zero energy to a maximum or end-point energy. From the end-point energy 
an expression (eq.. (l)) can he derived for a weight factor n(Ep), which, when 
plotted, graphically represents the beta spectrum. The beta spectrum can be 
divided into parts, each of which, as depicted in sketch (j), represents the per- 
centage of beta particles emitted in an energy segment and each of which is rep- 
resented by its midpoint energy Ep. Thus, in calculating n, the number of 
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particles reaching the collector, the calculation proceeds for each with 


N £ - Nn*(E p ) 


(where * refers to the normalized quantity) and then these partial results are 
summed over the appropriate normalizing factor. Hence, n takes the form 



— 


X™' n(E B ) 

p9 

I max 

/ r sin 0 cos 0 + t„(l - cos 0. )d0 

m max x 

AE 

E P 

+5 





^n(E p )AE 



E P 



(19) 
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In addition to the spectrum effect on the calculation, there exists another 
effect that must he considered at this point. As can he seen in sketch (k), 
beta particles emitted in the "backward" direction can reach the collector by- 



passing through the opposite side of the emitter (assuming no backscattering) . 
This effect can be handled in a rather simple way. These particles increase 
the number of particles in the forward direction by an amount just equal to 
an increase that would be caused by doubling the fuel thickness. Thus, in the 
calculations 


t f (calculations) = 2 times actual fuel-layer thickness (20) 

which will be used in all particle computations. The proof of this assertion 
rests on the fact that the two angles 0p made with the surface of the emitter 
by a particle described are equal for spherical and cylindrical configurations. 


Total Energy Reaching Collector 

If dn is the number of particles per second per unit emitter area reach- 
ing the col£ector from a volume element dv with an energy Ep^ c , the total 
energy per second per unit emitter area reaching the collector is 



(21) 


which becomes, after substituting for 


dn £ and summing over the energy spectrum. 




( 23 ) 


The energy E (r , 0) is determined, from 

V) c 

E = E . - (E + E + e$) 
p, c |3 f s 

where E^ and E g are the energies lost in the fuel and support layers, respec- 
tively. 

The quantities E-f and E s are found as follows (illustrated graphically): 
Sketch (l) represents a plot of range against energy for a beta particle in the 



Energy — ► 

fuel layer. For a travel distance of Ar in the fuel layer (from birth to 
arrival at the support layer), the beta particle will lose the energy 

Ej? ~ Eg — Eg (24:) 

On traveling through the support layer (see range energy curve in sketch (m)), 



the beta particle will lose the energy 

E = E 9 - E, 
s 2 o 


(25) 
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Hence ; on substituting in equations (24) and (25); equation (23) reduces to 

Ep, c (r,0) = E 5 - e$ (26) 

where E3 results from the preceding graphical calculations; and e$ is the 
energy the beta particles lose in the electric field. (This process can be re- 
peated until a set IL c ( r j_;®j[) i 8 generated.) 


Distribution of Particles at Collector 

From the integral expression for n e (eq. (7)) the following equation for 
An e can be written as 


N- 

An £ = — sin 0 cos 0 Ar (27) 

This equation represents the increment of particles leaving the emitter in a 
ringed segment of angle A9 and from a segment of fuel-layer thickness Ar. 

With the proper summation over the beta spectrum; equation (27) represents the 
distribution D(0, r) of particles leaving the emitter. Further; staying within 
the limits previously set forth (eq. 17)); equation (27) represents the distri- 
bution of Just those particles that reach the collector. 

To obtain the desired distribution, it is necessary to convert the distri- 
bution (eq. (27)) with its dependence on r and 0 to one with a dependence on 
Ep c , the energy of a particle at the collector, and y, the angle that the in- 
coming particle at the collector has with the normal. From the previous section 
it can be seen that for a given Ar^_ there exists a unique (Ep^ c ) . Hence for 
a ’ 1 

- ^,e > 1+1 < 28 > 


there exists a 


Ar = Ar^ - Ar.j_ +1 (29) 

Thus, by substituting AEp c for the corresponding Ar into equation (27), the 
distribution D(0, r) reeomes D(0,E c ), where E c is taken as the midpoint 

^e- 

From the equations of motion (eqs. (8) and (9)) set forth previously, an 
expression can be derived for 


r 


tan' 


-1 v t 


( 30 ) 


where Vj. and v r have previously been defined (sketch (f)). The expression 
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derived directly from equations (10a) and (lla) (or (lOb) and (lib)) with v r 
not set equal to zero is 


i , (M 2 

sin 2 y V R l/ sin 2 0 
or, nonrelativist ically, 

1 = /^2i\ 1 

. 2 \R, / . 2 a 

sin y \ 1/ sin 0 

1 2 

where IL (or 2 m o T o^ can determined from 

U 1 < 0r | m o v o } = E p,c + e$ < 32) 

For a given pair, 9 (midpoint of A9) and Ep c , there exists from equa- 
tion (3l) a y. Hence, the distribution D(0, E__ \ can be converted to 
D(y, Ep^ c ), where y is taken as the midpoint of Ay, by noting that 

^ = 0 i “ 9 i+l 

and 

* = r i ' r i + l 

Thus, in the preceding manner, the distribution of particles against angle and 
energy at the collector may be obtained. 

CYLINDRICAL CONFIGURATION 

Number of Particles Reaching Collector 

A spherical coordinate system for the cylindrical configuration is best set 
up as shown in sketch (n) (note that 0 is the complement of the usual 0). 

The volume element becomes 

dv ss r 2 cos 0 dr d0 dcp 
and the projected area becomes 

dS = dS cos 0 cos cp 



(U.. + U - S>e) 2 - U 2 
^ 1 o o 

( D 1 + U o> 2 - U o 
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Thus, following the argument set forth for the spherical geometry, the expres- 
sion for tog can he written as 

N 

dn' = — cos 9 cos m dS dr d8 dcp (33) 

e 4?t 



where n e is the number of particles reaching the collector per second per unit 
emitter area, and the limits are confined to one quadrant. These limits of in- 
tegration remain to be determined. 

For a beta particle in the electric field of the cylindrical cell, the 
equations of motion can be written as (see sketch (o)) 

Ui = Ug + eO (Conservation of energy) (35) 

1 = ^2^t 2 (Conservation of angular momentum) (36) 

p = p (Conservation of momentum) (37) 

z, 1 z, 2 

where the quantities in the equations (35) and (36) have been defined previously, 
and the quantity P* is the linear momentum in the z-direction. On substituting 
into the equations (35) to (37) they become 



16 



(39a) 



m o v l) min sln 9 


f- 


r \, min ' 


m v 
o z 


v t + T " 


(40a) 


where v r (radial Telocity at the collector) is taken as zero and t z i3 de- 
fined as the velocity of the beta particle in the z-direction at the collector. 
The equations can he solved to give 


(U. , + U - e$) 2 - U 2 

1, min o o 

l,mln * D o> 2 - D o 



2 2 2 
cos 9 sin cp + sin 9 


(41a) 


from which the minimum initial kinetic energy, can he found. (The non- 

relativistic formulations of the equations of motion (eqs. (38a) to (40a) are: 


(! m o T l,minj * | m o^t + T z) + e$ 
Vl,min cos 9 sln ^ R 1 = Vt 1 ^ 


(38b) 

(39h) 
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m Q v sin 9 = m 0 Y z 


(40b) 


On eliminating v^. and v z among them. 


m o v l, min) 
(I "4, min) 


e$ IR-i \ 2 ? 2 

= (_i) cos^0 sin cp + sin 0 

.R, 


(41b) 


is obtained. ) Continuing to follow the procedure previously set forth (for the 
spherical case), the limits r™ ax (0^cp), 0 Trinv (cp), and ( 0 = 0) can be ob- 

tained by replacing t g /cos 0 ^spherical case) by t g /(cos 9 cos cp) (cylindrical 
case). Sets of these limits can be produced by repetition of the generating 
process. 

It is again necessary to modify these limits by taking into account the geo- 
metric boundaries of the fuel layer. As a result, there exists a 0^ (cp) that 
satisfies the relation 


cos r (0 + ,cp)cos 


max' t ; 

t. 


<p 


i < 


~ r max (0 ^ )cos * 


0 t > 0 


0 t = 0 


and a cp. (0 = 0) that satisfies the relation 
*0 


008 h • r max (0/e t ) h > 0 


1 < 


- r (o,oT 
max' ’ 




= 0 


J 


Also, the values of 9 max and 0 max must not exceed rt/ 2. Therefore, 

^tj/cos 0 cos cp 


»«4 {/VV 


cos^0 cos cp dr 


/ ‘Pfc /*®ma x /* 1 max 

dcp / d 9 1 

*'04. •'0 


d9 / cos^0 cos cp dr 

'0 


(42) 


(43) 


/ ^max r ^max /* r max 

dcp / d0 / 


dcp / id f cos^d cos cp dr] 

'cp t Jo J 0 


(44) 
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which reduces to 




sin dcp + 


t 

Je\ 


dcp 



(?) 


L max 


cos 9 cos cp d 9 



(45) 


for beta particles of energy E p . Here again both forward directed and backward 
directed particles are accounted for by choosing t f as twice the actual fuel- 
layer thickness. 


Total Energy Reaching Collector 

From the discussion set forth for the spherical configuration, it can be. 
seen that the energy per second per unit emitter area reaching the collector for 
the cylindrical configuration is 



/ /* cp t f e t r 

n(E p ) / dcp / d 9 

VO ■'0 *T> 


.t^/cos 9 cos cp 


Ep c cos^0 cos cp dr 


/ 9 z»r 

max j max 

d V 

•'n 


dcp / d0 


cos 9 cos cp Ep^ c dr 


/ ^max /* 0 max /* r max \ 1—^ 

&9 I d el cos 2 0 cos cp Ep^ c drjAE jy , n(E B )AE 

n *4 •'o / / 


(46) 


The quantity Ep c (r, 9, cp) is determined in exactly the same manner as in the 
spherical case with the single change 


t t 

— - — (spherical case) -* ® (cylindrical case) 

cos 9 cos 9 cos cp 
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Distribution of Particles at Collector 
The quantity An , in the cylindrical configuration takes the form 


2 

An_ = — cos 0 cos cp Ar A0 Ap 

t jt 


( 47 ) 


deduced from the Integral expression (34). The distribution D(0,cp, r) obtained 
from An £ can be converted, as was done in the spherical case, into the form 
D ( 0, cp, Ep c ). The angle with the normal at which the beta particles enter the 
collector, y, is made up of the components 0 C and cp c defined as 


tan cp c = 


(48) 


and 


tan 0 C = 


v+ 


sin cp 


(49) 


By solving the equations of motion previously set forth for the cylindrical con- 
figuration (with v r not equal to zero), the following equations are obtained: 


X 

c 

R-, 


2 2 2 
sin <p c cos 0 sin cp 


(U. + U - «e) 2 - U 2 


(Ui + U 0 ) 2 - ^ 


- sin 01 


(50) 


1 _ 

1 

(U. + U - $e) 2 - U 2 
' 1 o o 

tan 2 0 

sin 2 0 

(»x + - U o 


- 1 


(51) 


or, nonrelativist ically. 


1 

W) 

, 

A 2 

m o v o - 

sln 2 cp c 

2n 

cos 0 sin 

2 

cp ' 

11 2 
\ 2 Vo 

1 

1 1 

(\ Vo - e$ 

tan 2 0 c 

sln 2 0 1 


1 m ,r 2 

2 ° O 


- sin 2 0| 


(50b) 


(51b) 
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( 52 ) 


where 


U-l (or 


1 2 n 
2 “W 


is determined from 


(or i m v(?) = E + e$ 

2 o o' p, c 

from which the angles cp c and 9 C can be deduced. Thus, with 


cos x = cos cp c cos 0 c 


(53) 


the distribution takes the form D(r, E ). 

v ' ’ p, c 


Efficiency of Cell 

For a cell operating at a voltage $, the electrical efficiency can be 
written as 


$en 
11 = t^E' 


(54) 


where n is the number of particles reaching the collector (per second), e the 
electron charge, tf the actual fuel- layer thickness, H the number of beta 
particles emitted per unit volume (per second), and E' their average energy. 
The average energy of the beta particles can be deduced from 


E' 



(55) 


where only the 2.98-Mev (end-point energy) beta spectrum is considered. The ex- 
pression for the efficiency (eq.. (54)) is valid for both spherical and cylin- 
drical configurations. 


Computer Calculations 

The computations were done with the IBM 7090 electronic computer following 
procedures (the relativistic equations and Katz-Penfold relations' were used) as 
outlined in the preceding sections. The programs (one for each of the two cell 
geometries) consisted of three major calculations: the number of particles 

reaching the collector (with the associated cell efficiency included in this 
section), the total energy of the particles arriving at the collector, and the 
distribution of particles in terms of energy and incoming angle at the collector. 

As an initial subordinate routine, the programs were set up to calculate 
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the beta spectrum weight factor n(Eo) (eq. (l).). It was found that answers of 
sufficient accuracy were obtained in the first two sections of the computer pro- 
grams whep the beta spectrum was divided into 15 segments, since the answers 
changed by less than 2 percent when twice (30) the number of segments were used. 
For the particle distribution calculations, the beta spectrum had to be divided 
into at least 30 segments to obtain smooth distributions. Integrations were 
performed on the computer by use of Simpson's rule. Since sets such as 
r max (0i) and Ep c (r^, 0^) could be obtained in the analysis, the use of numerical 
integration as performed on the computer was a convenient approach to the prob- 
lem. The number of steps (in the variables r, 0, and cp) taken in the integra- 
tion drastically affected the running time, and, therefore, a minimum number of 
steps, which allowed sufficient accuracy (less than 2-percent change on using 
half the step size), was used. The inverse of the Katz-Penfold relation 
(E = f(Rp)) was obtained through a polynomial curve fit in which terms beyond 
the order were neglected. 


DISCUSSION 

Table I is presented as a summary of the results of the first two sections 
of the computer programs. In the table a variety of cases set apart by their 
different radius ratios and layer thicknesses (expressed in cmXlO” 3 ) are pre- 
sented for both cell geometries. The cell parameters n, q, and E c are set 
down against operating voltage. The range of voltage is 0. 2 to 1.8 Mev taken in 
0. 2-Mev steps. The third section of the computer programs, which pertain to the 
particle distributions at the collector, is presented in graphical form. Plots 
representing pertinent features of the analysis drawn from the data of the 
tables are presented. In addition, curves are shown that compare the results of 
using different forms of the physical constraints. 

A discussion of cell design based on results of this analysis was presented 
in reference 6. 


Number of Particles Reaching Collector 

Curves of current (number of particles reaching the collector) against cell 
voltage are presented in figure 1 for different radius ratios and for both cell 
configurations. The significant feature of these curves is that an increase in 
the cell voltage corresponds to a decrease in the cell current. The curvature 
of these monotomic curves is entirely positive for all cylindrical cases, but a 
slight inverse curvature is detected in the higher radius-ratio cases of the 
spherical configuration (fig. 1(a)). This feature is probably the result of the 
form of the beta spectrum. An important aspect of these curves is that from 
them an operating point for the cell can be determined. In essence, an operat- 
ing point is reached when the cell and the attached device have matched (equal) 
voltages and current. 
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Efficiency 

Figure 2 represents curves of cell efficiency against cell voltage for 



various radius ratios and for both cell configurations. In the collection of 
monoenergetic particles (ref. 3), the maximum occurs at approximately half of 
the highest obtainable voltage (corresponding to the energy of the particle con- 
sidered). The effect of the beta spectrum changes this situation. The maximums 
in the plots occur between 0.4 and 0.9 megavolts, which is about half of the 
beta-particle average energy (1.31 Mev as determined from the computer). The 
long tails beyond 1.31 Mev are due to particles from the upper end of the beta 
spectrum reaching the collector. The weight per unit power output can be com- 
puted simply from the efficiency and weights of the various components. 


Energy Distribution 

Computations of the total energy of particles reaching the collector have 
been performed by using the analysis. Figure 3 presents typical curves. As the 
cell voltage increases, the number of particles reaching the collector rapidly 
diminishes, and thus the total energy also decreases. 

With knowledge of the cell efficiency (percent of energy initially available 
leaving the cell) and the percentage of the energy reaching the collector, the 
portion of the energy distribution absorbed In the components of the cell can be 
calculated. Plots showing this distribution of energy are presented for both 
configurations and typical thickness (fig. 4). It is assumed in these plots 
that the energy of the particles that does not reach the collector and is not 
part of the power output of the cell is completely absorbed (after what may be 
multiple passes) in the emitter foil (fuel and support layer). A knowledge of 
the distribution of energy in an operating cell Is important in determining com- 
ponent temperatures of the cell (see ref. 6). 


Particle Distributions at Collector 

Since the impingement of high-energy electrons on the collector can cause 
secondary effects, the study of the distribution of particles at this component 
is an important consideration. Particle distributions E(y, E^ c ) 'that represent 
the number of particles arriving at the collector as functions of particle energy 
and arrival angle have been evaluated for a single case for both configurations. 
From this information, two types of curves have been plotted: (l) The number of 

particles in an increment of energy and angle against arrival angle (midpoint) 
at a fixed particle energy increment was drawn and (2), the same dependent vari- 
able against particle energy (midpoint) for a fixed increment of arrival angle 
was also drawn. The increments of energy and angle taken were 0.1 Mev and 5°, 
respectively. The operating voltage is 0.6 Mev for both cases considered. The 
plots are proportional to § 2 n/dEp Q hr (the second partial of n with respect 
to particle energy and arrival angle). These two types of curves are shown in 
figure 5. 

The calculated points did not represent perfectly smooth curves. Smooth- 
ness of the distributions is a function of step size in the computer program. 
Running time is also a function of step size, but, because smoothness and a short 
running time are not compatible, a compromise was made. 
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Some of the features of figure 5 may now he noted. It is observed for the 
spherical case (figs. 5(a) to (e)) that the particles arriving with a high 
energy cannot have an arrival angle greater than some "cutoff" angle, which dif- 
fers for different particle arrival energies. It is further noted for the 
spherical case that particles arriving in a fixed increment of arrival angle 
tend to favor a small range of energies (see figs. 5(f) to (h)). In the cylin- 
drical case (figs. 5(i) to (o)), this bunching is not observed in either type of 
plot. Instead, the curves spread out over the entire range of the variables con- 
sidered. The nature of the spreading is a result of the geometry, but the fact 
that spreading does exist is due mainly to the beta spectrum and slightly to the 
different origins of the beta particles in the fuel layer. 


Effects of Changes in Physical Relations 

Figure 6 shows an interesting comparison of cell efficiency plots for the 
cerium-praseodymium 2. 98-Mev beta particle and a hypothetical 2.18-Mev beta 
particle from an allowed nuclear transition. Changing the beta spectrum not 
only affects the spread of the curve but also affects the value of the maximum. 

Figure 7 compares the results of using the relativistic and nonrelativlstlc 
equations in a typical efficiency calculation. As would be expected, the effect 
of the relativistic equations becomes pronounced at high energies, but the dif- 
ference between the two efficiency calculations is small at optimum cell voltage. 

Figure 8 compares the Katz-Penfold range-energy relation for cerium and 
molybdenum (a possible support-layer material) with a linear approximation. 
Results of using these two different forms of the range-energy relation are com- 
pared in figure 9 for a typical efficiency calculation. The difference is seen 
to be quite small. 


CONCLUDING REMARKS 

A mathematical analysis of the beta radioisotope cell has been developed 
for calculations with an electronic computer. The beta radioisotope considered 
was cerium 144, which in the decay chain yields a beta-particle spectrum with a 
2. 98-Mev end-point energy. Results have been discussed and graphically illus- 
trated. 

From the results of the calculations, certain conclusions can be drawn. 

For a given radius ratio and component thickness, the maximum cell efficiency 
is significantly greater in the spherical configuration than in the correspond- 
ing cylindrical configuration. It is also observed that cell efficiency in- 
creases with radius ratio and that the optimum voltage (corresponding to maximum 
efficiency) also increases. In particle distributions at the collector it is 
found that cutoff angles exist in the spherical configuration where none are ob- 
served for the cylindrical configuration. Furthermore, it is found that for the 
lower arrival energies more particles are impinging on the collector with a 
grazing angle. 
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Changing the beta spectrum drastically affects the cell efficiency. The 
comparison of the relativistic equations -with the nonrelativist ic equations shows 
a small hut increasing effect at high cell potentials. Finally, the use of a 
linear approximation in the range-energy relation has a small effect on the cal- 
culations. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 25, 1963 
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APPENDIX - SYMBOLS 


A 

C 

D ^ E p,c> 

E 

E 

c 

E f 

E o 
E p> c 
E s 
Ep 

e 

m o 

N 

n 


n(Ep) 

n*(E p ) 





r 


dr 


atomic mass number 
speed of light 
distribution at collector 
energy 

total energy arriving at collector 
energy lost in fuel layer 
end-point energy 

energy of beta particle at collector 
energy lost in support layer 
energy of emitted beta particle 
electronic charge 
mass of electron 

number of disintegrations per unit volume per second 

number of particles per unit emitter area per second reaching col- 
lector 

relative number of particles in beta spectrum 

normalized n(Ep) 

linear tangential momentum 

linear momentum in z-direction 

practical range; effective straight-line distance in -which an elec- 
tron loses its entire initial energy 

inner radius 

outer radius 

radial distance 

differential element of radial distance 
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r B maximum radial distance 
dS differential element of surface area 

tj, fuel-layer thickness 

t s support-layer thickness 

U kinetic energy of beta particle 

U 0 rest mass energy 

kinetic energy of beta particle emitter 

U 2 kinetic energy of beta particle collector 

V operating voltage, Mev 

v velocity 

v_ velocity of beta particle in z-dlrection at collector 

Z atomic number 

T arrival angle (normal) at collector 

e defined in eq_. (l) 

T) cell efficiency 

8 polar angle 

d Q 0-component of arrival angle 

maximum polar angle 

dv differential element of volume 

$ cell voltage 

qp azimuthal angle 

cp c cp-component of arrival angle 

cp maximum azimuthal angle 

ihhx 

Superscript 


T 


times a unit area 



Subscripts: 


max 

min 

P 

e 

1 

2 


maximum 

minimum 

projected area 

in an increment of energy 

initial (emitter) 

final (collector) 
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TABLE I 


COMPUTER OUTPUT FOR CERIUM 144 CELL WITH MOLYBDENUM SUPPORT LAYER 


(a) Spherical configuration 




Fuel 

-layer thickness, t f , cm 



1. 27X10" 3 

2. 5 4X10“ 3 



Support-layer thickness, 

b s , cm 



2. 5 4X10“ 3 

2. 5 4x10" 3 




Radius ratio, Rg/R 

1 



2.0 

2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency, 

energy 


n, 


arriving 


n, 


arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




E c> + 
percent 




E c> 

percent 

0.2 

83.8 

12.8 

66.8 

0.2 

82.1 

12.5 

64.4 

. 4 

74.9 

22.9 

53.9 

. 4 

73.2 

22.3 

51.8 

.6 

62.6 

28.6 

41.5 

.6 

61.2 

28.0 

39.8 

.8 

48.6 

29. 7 

30.1 

.8 

47.1 

28.8 

28.9 

1.0 

33.9 

25.9 

20.2 

1.0 

32.4 

24.8 

19.1 

1.2 

21.5 

19.7 

12.0 

1.2 

20.3 

18.6 

11.2 

1.4 

11.4 

12.2 

5.85 

1.4 

11.1 

11.9 

5.53 

1. 6 

5.80 

7.09 

2.50 

1.6 

5.81 

7.10 

2. 43 

1.8 

2.72 

3.74 

.96 

1.8 

2.73 

3. 75 

.93 

Fuel-layer thickness, tf, cm 

5.08X10 -3 

5. 08xl0“ 3 



Support-layer thickness, 

t s , cm 



2. 5 4X10“ 3 

2.54X10 -3 

Radius ratio, R 2 /R 1 


1 . 

L 


1.3 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles, 

ficiency, 

energy 

voltage, 

particles. 

ficiency, 

energy 


n. 

T\, 

arriving 


n. 

ib 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




E c> 

percent 




E c> 

percent 

0.2 

58.6 

8.95 

51.8 

0.2 

73.5 

11.2 

59.2 

.4 

35.0 

10.7 

28.9 

.4 

53. 6 

16.3 

43.1 

.6 

22.0 

10.1 

16.5 

.6 

34.2 

15.6 

26.1 

.8 

14.0 

8.55 

9. 60 

.8 

20.8 

12.7 

14.2 

1.0 

8.64 

6.59 

5.30 

1.0 

12.7 

9.70 

7.85 

1.2 

5.07 

4.65 

2.76 

1.2 

7.40 

6.77 

4.07 

1.4 

2.66 

2.84 

1.25 

1.4 

4. 02 

4.29 

1.95 

1.6 

1.48 

1.80 

.61 

1.6 

2.12 

2.60 

.86 

1.8 

.61 

.84 

.20 

1.8 

.94 

1.29 

.23 




TABLE I. - Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 
WITH MOLYBDENUM SUPPORT LAYER 
(a) Continued. Spherical configuration 



Fuel- layer thickness, tf, cm 


5.08x10-3 


5.08x10-3 



Support- layer thickness, t , cm 


2.54x10-3 2. 54x10“ 3 


Radius ratio, R 2 /R 1 


Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles. 

ficiency, 

energy 

voltage, 

particles, 

ficiency. 

energy 

®, 

n, 

T)> 

arriving 

' ®, 

n, 

lb 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




®c> 




®c> 




percent 




percent 

0.2 

75.5 

11.5 

59. 6 

0.2 

78. 3 

11.9 

60.0 

.4 

62.0 

18.9 

46.4 

.4 

66.9 

20.4 

47.6 

.6 

45.5 

20.8 

32.9 

.6 

54. 2 

24.8 

36.0 

.8 

29.3 

17.8 

20.2 

.8 

39.2 

23.9 

25.0 

1.0 

17. 7 

13.5 

11.0 

1.0 

25. 6 

19.6 

15.5 

1.2 

10.1 

9.22 

5.55 

1.2 

15.1 

13.8 

8. 40 

1. 4 

5.41 

5.78 

2.62 

1.4 

8.20 

8. 77 

3. 95 

1. 6 

2.83 

3. 45 

1.13 

1. 6 

4. 18 

5.10 

1. 70 

1.8 

1.33 

1.83 

.46 

1.8 

1. 92 

2.65 

.64 


Fuel- layer thickness, t f , 


5. 08X10“ 3 


5.08x10" 


Support- layer thickness, t s , cm 
2. 54x10“ 3 2. 54x10“ 3 

Radius ratio, R 2 /R]_ 


Cell Number of Cell ef- 
voltage, particles, ficiency, 
n, n> 

mega- percent percent 
volts 


Total 
energy 
arriving 
at col- 
lector, 

E c , 

percent 


Cell Number of Cell ef- 
voltage, particles, ficiency, 

mega- percent percent 
volts 














TABLE I. 


Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 


WITH MOLYBDENUM SUPPORT LAYER 
(a) Continued. Spherical configuration 





Fuel- 

layer thickness, tf, cm 




5. 

08X10- 3 


5.08x10-3 




Support-layer thickness, 

t s , cm 




2. 

5 4X10" 3 


1. 27X10" 3 

Radius ratio, Rg/R^ 

10.0 

1.5 

Cell 

Number 

of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency, 

energy 

®, 

n. 


ib 

arriving 

®> 

>b 

ib 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 




lector, 

volts 



lector, 





percent 




percent 

0.2 

79.5 


12.1 

60.0 

0.2 

78.9 

12.0 

63.7 

. 4 

72.7 


22.2 

48.2 

.4 

65.1 

19.7 

50.0 

.6 

65.5 


30.0 

37.5 

.6 

47.5 

21.8 

35. 6 

.8 

56.5 


34.5 

28.1 

.8 

30.9 

18.8 

21.9 

1.0 

46.6 


35.6 

20.0 

1.0 

18.5 

14.1 

11.9 

1.2 

36.3 


33.2 

13.8 

1.2 

10.5 

9. 62 

5.89 

1.4 

26.4 


28.2 

8.81 

1.4 

5.88 

6.28 

2.89 

1.6 

18.0 


22.0 

5.09 

1.6 

3.06 

3.74 

1.26 

1.8 

11.4 


15.6 

2.54 

1.8 

1. 45 

1. 99 

.49 




Fuel- 

•layer thickness, t-f, cm 




5. 

08X10" 3 


5.08x10-3 




Support- layer thickness, 

; s , cm 




1. 

27x10-3 


1. 27x10-3 





Radius ratio, Rg/R 

1 





2. 

0 


4.0 

Cell 

Number 

of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency, 

energy 

®, 



Tb 

arriving 

®, 

n, 

ib 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 




lector, 

volts 



lector, 





®c> 

percent 




percent 

0.2 

80.9 


12.4 

64.0 

0.2 

82.2 

12.5 

64.0 

.4 

72.0 


22.0 

51.5 

.4 

75.3 

23.0 

51.7 

.6 

59.9 


27.4 

39.6 

.6 

66.7 

30.5 

40.4 

.8 

45.9 


28.0 

28.7 

.8 

57.0 

34.8 

30. 4 

1.0 

31.6 


24.1 

18.9 

1.0 

46. 4 

35.4 

22.0 

1.2 

19.9 


18.2 

11.0 

1.2 

34.5 

31.6 

15.1 

1.4 

11.3 


12.0 

5.58 

1.4 

24.2 

25.8 

9.51 

1.6 

5.67 

6.92 

2.37 

1.6 

15.5 

18.9 

5.38 

1.8 

2.70 

3.71 

.92 

1.8 

8.88 

12.2 

2.62 




TABLE I. - Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 
WITH MOLYBDENUM SUPPORT LAYER 
(a) Continued. Spherical configuration 




Fuel- layer thickness, 1 

f, cm 



5. 08x10“ 3 

5. 08X10“ 3 



Support-layer thickness, 

t s , cm 



5. 08X10“ 3 

5. 08xl0“ 3 

Radius ratio, Ro/Ri 


1 . 

5 


2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles. 

f iciency. 

energy 

voltage, 

particles. 

f iciency, 

energy 


n, 


arriving 

®, 

n, 


arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent. 

at col- 

volts 



lector, 

volts 



lector, 








®c> 




percent 




percent 

0.2 

33.3 

5.08 

20.5 

0.2 

72.1 

11.0 

53.0 

.4 

25.8 

7.88 

14.7 

.4 

63.3 

19.3 

42.8 

.6 

18.1 

8.30 

9.52 

.6 

51.7 

23. 7 

31.7 

.8 

11.5 

7.04 

5.47 

.8 

38.9 

23.7 

22.5 

1.0 

6.78 

5.17 

2.81 

1.0 

26.4 

20. 4 

14.6 

1.2 

3.74 

3.43 

1.30 

1.2 

16.2 

14.8 

8.28 

1.4 

1.94 

2.07 

.56 

1. 4 

9.03 

9. 65 

4. 12 

1.6 

.90 

1.10 

.21 

1. 6 

4. 66 

5. 70 

1. 78 

1.8 

.38 

.53 

.06 

1.8 

2.10 

2.88 

.65 



Fuel- layer thickness, t 

f , cm 



5. 08x10“ 3 

10. 16X10“ 3 



Support- layer thickness, 

t s , cm 



10. 16X10“ 3 

2. 54x10“ 3 

Radius ratio, R 2 /R 1 

2.0 

2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

f iciency. 

energy 

voltage, 

particles, 

ficiency, 

energy 




arriving 


n, 


arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




*0 








percent 




percent 

0.2 

62.8 

9.59 

42. 5 

0.2 

70.3 

10.7 

52. 4 

. 4 

53.9 

16.4 

33.0 

. 4 

60.8 

18.5 

41. 4 

.6 

43.8 

19.8 

24.3 

. 6 

50.0 

22.9 

31.2 

.8 

31.4 

19.2 

16.7 

.8 

37. 6 

29.9 

22.0 

1.0 

20.8 

15.8 

10. 4 

1.0 

25. 7 

19. 6 

14.2 

1.2 

12.5 

11.5 

5.82 

1.2 

15.8 

14.5 

8.06 

1.4 

6.89 

7.36 

2.79 

1.4 

8.82 

9. 43 

3. 96 

1.6 

3.43 

4.19 

1. 15 

1.6 

4.57 

5.58 

1. 72 

1.8 

1.52 

2.09 

.40 

1.8 

2.12 

2.92 

.64 
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TABLE I. 


Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 


WITH MOLYBDENUM SUPPORT LAYER 
(a) Continued. Spherical configuration 


Fuel-layer thickness, tf, cm 

12. 7x10" 3 

20.24X10- 3 



Support-layer thickness. 

t s , cm 



5. 08X10“ 3 

2. 54xl0“ 3 

Radius ratio, Rg/Ri 


2. 

0 


2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency, 

energy 


n, 

»b 

arriving 



ib 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




E c , 




®e. 




percent 




percent 

0.2 

62.2 

9.50 

43.9 

0.2 

58.1 

8.87 

41. 4 

.4 

53. 6 

16.8 

34.2 

.4 

50.0 

15.3 

32.1 

. 6 

43.3 

19.8 

25.4 

.6 

40.3 

18.5 

23. 7 

.8 

31.9 

19.5 

17.5 

.8 

30.0 

18.3 

16.3 

1.0 

21.3 

16.2 

11.0 

1.0 

20.0 

15.2 

10.2 

1.2 

13.0 

11.9 

6.10 

1.2 

12.2 

11.2 

5.64 

1.4 

7.22 

7.72 

2.96 

1.4 

6.75 

7.21 

2.72 

1.6 

3.64 

4. 45 

1.23 

1.6 

3.42 

4.17 

1.11 

1.8 

1. 64 

2.25 

.44 

1.8 

1. 55 

2.14 

.39 

Fuel- layer thickness, t f , cm 

25.4X10- 3 

40. 64x10-3 



Support- layer thickness, 

; s , cm 



5.08X10- 3 

2.54x10-3 

Radius ratio, Rg/Ri 

2.0 

2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency, 

energy 


n, 


arriving 

*> 

n, 

ib 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector. 












percent 




percent 

0.2 

49.9 

7. 62 

33.2 

0.2 

42.3 

6.46 

27. 6 

.4 

42.2 

13.0 

25.2 

.4 

35.5 

10.8 

20.8 

.6 

33.5 

15.3 

18.2 

.6 

28.0 

12.8 

14.8 

.8 

24.3 

14.8 

12.3 

.8 

20. 3 

12. 4 

9.91 

1.0 

16.0 

12.2 

7. 47 

1.0 

13. 4 

10.2 

5.98 

1.2 

9.50 

8. 70 

3. 97 

1.2 

7.99 

7.32 

3.18 

1. 4 

5.12 

5.54 

1.87 

1.4 

4. 32 

4. 61 

1. 46 

1.6 

2.58 

3.15 

.74 

1.6 

2.13 

2.60 

.57 

1.8 

1.13 

1.55 

.25 

1.8 

.93 

1.28 

.20 




TABLE X. - Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 
WITH MOLYBDENUM SUPPORT LAYER 
(a) Concluded. Spherical configuration 


101. 6X10" 


5.08x10" 


Fuel- layer thickness, tj, cm 

0.12 7X10“ 3 

Support- layer thickness, t s , cm 


2.54x10" 


Radius ratio, Rg/Ri 



1.5 




2.0 



Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles. 

ficiency. 

energy 

voltage, 

particles, 

ficiency, 

energy 

®, 

n, 

ib 

arriving 

a>> 

n, 

lb 

arriving 

mega- 

volts 

percent 

percent 

at col- 
lector, 

percent 

mega- 

volts 

percent 

percent 

at col- 
lector, 

percent 

0.2 

18. 6 

2.83 

10.5 

0.2 

84.3 

12.7 

69.1 

. 4 

14.2 

4. 34 

7.39 

. 4 

75.7 

23.1 

55.8 

.6 

9.84 

4.50 

4.70 

.6 

63.5 

29.1 

43.1 

.8 

6.17 

3.77 

2.66 

.8 

48.7 

29.7 

31.1 

1.0 

3.57 

2.72 

1. 33 

1.0 

34.3 

26.2 

20.1 

1.2 

1.93 

1.76 

.68 

1.2 

21.5 

19.7 

12.4 

1.4 

.97 

1.04 

.26 

1.4 

11.5 

12.3 

5.99 

1.6 

.45 

.55 

.09 

1.6 

5.80 

7.09 

2.57 

1.8 

.19 

.26 

.03 

1.8 

2.72 

3.74 

.99 
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TABLE I. - Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 


WITH MOLYBDENUM SUPPORT LAYER 
(t>) Cylindrical configuration 


Fuel- layer thickness, tf, cm 

1.27X10" 3 

2.54X10" 3 

Support- layer thickness, t g , cm 

2. 54x10" 3 

2. 54x10" 3 

Radius ratio, R 2 /R 1 

2.0 

2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles, 

ficiency, 

energy 

voltage, 

particles. 

ficiency, 

energy 


b 


arriving 

V 

*b 

b 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector. 

volts 



lector, 




®c, 




E c , 




percent 




percent 

0.2 

68.1 

10.4 

58.4 

0.2 

67.3 

10.2 

57.0 

.4 

51.6 

15. 7 

40.6 

. 4 

50.8 

15.5 

39.8 

. 6 

37.1 

16.9 

26. 9 

• 6 

36.8 

16.8 

26. 6 

.8 

25.8 

15.7 

16.6 

.8 

25.0 

15.2 

16.6 

1.0 

16.8 

12.8 

9.59 

1.0 

15.9 

12.1 

9.53 

1.2 

9.83 

8.97 

5.10 

1.2 

9.22 

8.42 

5.04 

1. 4 

5.01 

5. 33 

2. 42 

1.4 

4.83 

5.15 

2.36 

1.6 

2.32 

2.83 

1.00 

1.6 

2.35 

2.86 

.98 

1.8 

1.02 

1.40 

.37 

1.8 

1.04 

1.42 

.37 

Fuel-layer thickness, t f , cm 

5. 08X10" 3 

5. 08x10" 3 



Support-layer thickness, 

t s , cm 



2. 54X10" 3 

2.54X10" 3 




Radius ratio, R 2 /R 

1 




1 . 

1 



1.3 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency, 

energy 


n. 

ib 

arriving 

*, 

n, 

h 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




®c> 




Ec> 




percent 




percent 

0.2 

51.4 

7.82 

44.8 

0.2 

59. 6 

9.07 

50.0 

.4 

30.5 

9.28 

24.8 

.4 

39. 4 

12.0 

31. 6 

. 6 

19.2 

8.77 

14.2 

.6 

24.2 

11.1 

18.2 

.8 

11.9 

7.29 

8.04 

.8 

14. 7 

8.93 

9.92 

1.0 

7.24 

5.50 

4. 39 

1.0 

8. 77 

6.67 

5.31 

1.2 

4.18 

3.81 

2.26 

1.2 

4. 95 

4. 52 

2. 71 

1.4 

2.23 

2.38 

1.06 

1.4 

2.86 

3.04 

1.32 

1.6 

1. 22 

1.48 

.48 

1.6 

1.47 

1.80 

.57 

1.8 

.50 

.67 

.17 

1.8 

. 64 

.88 

.21 




T ABU S I. - Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 
WITH MOLYBDENUM SUPPORT LAYER 
(b) Continued. Cylindrical configuration 


Fuel-layer thickness, tf, cm 


5.08x10-3 5.08x10-3 


Support-layer thickness, t s , cm 


2 . 5 4x10“ ^ 2 . 54x10" 3 


Radius ratio, Ro/Ri 


1.8 


Fuel-layer thickness, tf, cm 


5.08x10" 


2.54x10" 


5.08x10" 


Support- layer thickness, t s , cm 


2.54x10" 


Radius ratio, R 2 /R 1 


®C; 

percent 



4.0 



Cell 

Number of 

Cell ef- 

Total 

■voltage, 

particles, 

ficiency, 

energy 




arriving 

mega- 

volts 

percent 

percent 

at col- 
lector, 
®c> 

percent 

0.2 

66.9 

10.2 

53.5 

.4 

52.4 

16.0 

38.0 

.6 

40.0 

18.2 

26.1 

.8 

29.2 

17.7 

17.0 

1.0 

20.2 

15.4 

10.6 


13.6 

12.3 

6.20 

wSm 

8.52 

9.07 

3.34 


4. 72 

5. 75 

1.58 

1.8 

2.12 

2.99 

.63 

























TABLE I 


Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 


WITH MOLYBDENUM SUPPORT LAYER 

(t>) Continued. Cylindrical configuration 




Fuel- layer thickness, t 

f, cm 



5. 08X10" 3 

5. 08X10" 5 



Support-layer thickness, 

t s , cm 



2.54xi0" 3 

1.2 7X10" 3 

Radius ratio, R?/®! 


10. 

0 


1.5 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage. 

particles, 

ficiency, 

energy 



lb 

arriving 

■®, 

n, 


arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector. 

volts 



lector. 




percent 




E C 

percent 

0.2 

67.4 

10.3 

53.8 

0.2 

64. 7 

9.85 

55.0 

. 4 

53.1 

16.2 

38.3 

.4 

46.2 

14.1 

37.0 

.6 

40.9 

18.7 

26.4 

. 6 

30.5 

13.9 

23.2 

.8' 

30. 4 

18.5 

17.4 

.8 

19.0 

11.5 

13. 7 

1.0 

21.8 

16.5 

10.9 

1.0 

11.0 

8.38 

6.87 

1.2 

14.7 

13.4 

6.53 

1.2 

6.25 

5.70 

3. 46 

1.4 

9. 76 

10.4 

3. 62 

1.4 

3. 42 

3. 63 

1. 65 

1.6 

5. 95 

7.24 

1.79 

1.6 

1. 72 

2.09 

.70 

1.8 

3.12 

4.38 

• 75 

1.8 

.79 

1.08 

.27 

Fuel- layer thickness, t^, cm 

5. 08x10" 3 

5.08X10" 3 



Support-layer thickness. 

; s , cm 



1.2 7X10" 3 

1. 2 7X10" 3 

Radius ratio, 1 * 2 /®! 


2. 

0 


4.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency. 

energy 

voltage, 

particles, 

ficiency, 

energy 


n, 

% 

arriving 

*, 


ib 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




percent 




Ec 

percent 

0.2 

67.2 

10.2 

56.3 

0.2 

69.1 

10.5 

57.2 

.4 

51.0 

15.5 

39.5 

.4 

54.3 

16.5 

40.9 

.6 

37.0 

16.9 

26.4 

.6 

41. 6 

19.0 

28.3 

.8 

25.4 

15.5 

16.7 

.8 

30.8 

18.7 

18.8 

1.0 

15.9 

12.1 

9.87 

1.0 

21.7 

16.5 

11.8 

1.2 

9.31 

8.50 

5.30 

1.2 

14.3 

13.1 

7.08 

1.4 

5.01 

5.33 

2.46 

1.4 

9.03 

9.62 

3. 91 

1.6 

2.42 

2.95 

.99 

1.6 

5.18 

6.30 

1. 91 

1.8 

1.11 

1.52 

.37 

1.8 

2.57 

3. 53 

. 77 




TABLE I 


Continued. COMPUTER OUTPUT FOR CERIUM 144 CELL 


WITH MOLYBDENUM SUPPORT LAYER 
(t>) Continued. Cylindrical configuration 




Fuel- 

-layer thickness, t f , cm 



5.08x10-3 

5.08 X 10 -3 

Support- layer thickness, t s , cm 

5.08 X 10 -3 

5.08 X 10" 3 

Radius ratio, Rg/ E l 


1 . 

5 


2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles, 

ficiency. 

energy 

voltage, 

particles. 

ficiency, 

energy 

®, 

n. 


arriving 


n, 

ib 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




®c> 




®c> 




percent 




percent 

0.2 

29.3 

4. 46 

17.7 

0.2 

60.7 

9.24 

47.0 

.4 

19.8 

6.02 

11.2 

.4 

45. 6 

13.9 

32.4 

.6 

12.5 

5.70 

6.52 

.6 

32.4 

14.7 

21.2 

.8 

7.35 

4. 47 

3. 45 

.8 

21.4 

13.1 

13.0 

1.0 

4.07 

3.10 

1.68 

1.0 

13.2 

10.0 

7.37 

1.2 

2.13 

1.95 

.76 

1.2 

7.50 

6.85 

3.83 

1.4 

1.03 

1.10 

.30 

1.4 

4.04 

4.30 

1.81 

1.6 

.45 

.55 

.11 

1.6 

1.98 

2.41 

.75 

1.8 

.17 

.24 

.03 

1.8 

.82 

1.13 

.25 



Fuel 

-layer thickness, tf, cm 



5.08 X 10-3 

10.16 x 10-3 

Support- layer thickness, t s , cm 

10. 16 X 10-3 

2. 54 x 10-3 

Radius ratio, R 2 /R 1 


2. 

0 


2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage. 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency. 

energy 



Tb 

arriving 


n, 

Tb 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector, 




E c> 




E c > 




percent 




percent 

0.2 

52.8 

8.03 

38.2 

0.2 

60. 4 

9.18 

46. 4 

.4 

38.8 

11.7 

25.9 

,4 

45.1 

13.7 

32.0 

.6 

26.8 

12.2 

16.6 

.6 

32.2 

14.6 

20.9 

.8 

17. 4 

10. 6 

10.0 

.8 

21.7 

13.2 

12.9 

1.0 

10.7 

8.14 

5.59 

1.0 

13.4 

10.2 

7. 31 

1.2 

5.87 

5.36 

2. 77 

1.2 

7. 73 

7.06 

3.81 

1.4 

3.13 

3.33 

1.27 

1.4 

4.18 

4. 45 

1.79 

1.6 

1.52 

1.85 

.51 

1.6 

2.06 

2.50 

.75 

1.8 

. 62 

.84 

.16 

1.8 

.91 

1.25 

.27 




TABUS I. - Concluded. COMPUTER OUTPUT FOR CERIUM 144 CELL 
WITH MOLYBDENUM SUPPORT LAYER 
(b) Concluded. Cylindrical configuration 


Fuel- layer thickness, t f , cm 

12. 7X10“ 3 

20. 2 4X10“ 3 



Support- layer thickness. 

t s , cm 



5. 08x10“ 3 

2. 54x10“ 3 

Radius ratio, Rs/^l 

2.0 

2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency, 

energy 



ib 

arriving 


n. 

n> 

arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector, 

volts 



lector. 




=0* 




E c> 




percent 




percent 

0.2 

54.1 

8.23 

39.2 

0.2 

51.5 

7.84 

36.8 

.4 

40.1 

12.2 

26.8 

.4 

37. 9 

11.5 

25.1 

.6 

28.2 

12.8 

17.4 

.6 

26. 6 

12.1 

16.3 

.8 

18.5 

11.2 

10.5 

.8 

17.5 

10.6 

9.75 

1.0 

11.3 

8.63 

5.90 

1.0 

10. 7 

8.15 

5.45 

1.2 

6.51 

5.94 

3.00 

1.2 

6.10 

5.57 

2.77 

1.4 

3.41 

3. 63 

1.40 

1.4 

3.18 

3.39 

1.24 

1.6 

1. 67 

2.03 

.56 

1.6 

1.51 

1.84 

.49 

1.8 

.69 

.95 

.18 

1.8 

.64 

.87 

.16 

Fuel- layer thickness, tf, cm 

25. 4x10“ 3 

40. 64xl0“ 3 



Support-layer thickness, 

t s , cm 



5. 08X10“ 3 

2. 54x10“ 3 

Radius ratio, R 2 /R;[_ 

2.0 

2.0 

Cell 

Number of 

Cell ef- 

Total 

Cell 

Number of 

Cell ef- 

Total 

voltage, 

particles, 

ficiency, 

energy 

voltage, 

particles, 

ficiency. 

energy 

*, 

n. 


arriving 




arriving 

mega- 

percent 

percent 

at col- 

mega- 

percent 

percent 

at col- 

volts 



lector. 

volts 



lector. 




E c> 




E c; 




percent 




percent 

0.2 

44.4 

6.75 

29.7 

0.2 

37.9 

5.77 

24.3 

. 4 

32.3 

9.83 

20.0 

.4 

27.3 

8.30 

16.2 

.6 

22.2 

10.2 

12.7 

.6 

18.6 

8.52 

10.2 

.8 

14.3 

8.72 

7.51 

.8 

11.9 

7.26 

5.94 

1.0 

8.60 

6.54 

4.08 

1.0 

7.07 

5.37 

3.17 

1.2 

4. 77 

4.36 

2.00 

1.2 

3.85 

3. 52 

1.52 

1.4 

2.45 

2.61 

.89 

1.4 

1.93 

2.06 

.65 

1.6 

1.14 

1.38 

.33 

1.6 

.87 

1.06 

.24 

1.8 

.45 

.62 

.10 

1.8 

.34 

.47 

.07 
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Cell voltage, megavolts 
(a) Spherical configuration. 


Figure 1. - Humber of particles reaching collector as a function of cell voltage. Fuel-layer thickness, 35.0 milligrams 
per square centimeter; support-layer thickness, 25.9 milligrams per square centimeter. 





Cell efficiency, percent 



mm 



Cell efficiency, percent 




Cell voltage, megavolts 
(a) Spherical configuration. 


Figure 3. - Total energy reaching collector as a function of cell voltage. Fuel- layer thickness, 35.0 milligrams per 
square centimeter: support-layer thickness, 25.9 milligrams per square centimeter; radius ratio, 2.0. 



.2 


.4 


.6 


Cell voltage, megavolts 
(h) Cylindrical configuration. 

Figure 3. - Concluded. Total energy reaching collector as a function of cell voltage. Fuel-layer thickness, 35.0 mill! 
grains per square centimeter; support-layer thickness, 25.9 milligrams per square centimeter; radius ratio, 2.0. 


Cell voltage, megavolts 
(a) Spherical configuration. 


Figure 4. - Energy distribution in various cell components. Fuel-layer thickness, 70.0 milligrams per square 
centimeter; support-layer thickness, 25.9 milligrams per square centimeter; radius ratio, 2.0. 








(b ) Cylindrical configuration. 

Figure 4. - Concluded. Energy distribution in various cell components. Fuel-layer thickness, 70.0 milligrams 
per square centimeter; support-layer thickness, 25.9 milligrams per square centimeter; radius ratio, 2.0. 










Arrival angle at collector, y> deg 

(l) Cylindrical configuration; energy, 1.6 to 1.7 Mev. 

Figure 5. - Continued. Particle distribution at collector. Fuel-layer thickness, 55.0 milligrams 
per square centimeter; support-layer thickness, 25.9 milligrams per square centimeter; radius 
ratio, 2.0; cell voltage, 0.6 megavolts. 






Number of particles, percent 



(m) Cylindrical configuration; arrival angle 20° to 25°. 




0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 

Particle energy, Mev 


(o) Cylindrical configuration; arrival angle, 60° to 65°. 

Figure 5. - Concluded. Particle distribution at collector. 
Fuel-layer thickness, 35.0 milligrams per square centimeter; 
support-layer thickness, 25.9 milligrams per square centi- 
meter, radius ratio, 2.0; cell voltage, 0.6 megavolts. 
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Cell efficiency, percent 


Nonrelat ivist ic 


Relativistic 


Cell voltage, megavolts 

Figure 7. - Comparison of relativistic against nonrelat ivist ic formulations in cell efficiency calculation for 
spherical configuration. Fuel-layer thickness, 35.0 milligrams per square centimeter; support-layer thickness 
12.9 milligrams per square centimeter; radius ratio, 4.0. 












NASA-Langley, 1963 E-2108 


U1 

Ol 



Figure 9. - Comparison of different range-energy relations In efficiency calculations for spherical configura- 
tion. Fuel-layer thickness, 35. 0 milligrams per square centimeter; support -layer thickness, 12. 9 milligrams 
per square centimeter; radius ratio, 4.0. 



